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Abstract 

Th e empi rical copula process plays a central role for statistical inference on copulas. Recently, 



Seeersl (|201lh investigated the asymptotic behavior of this process under non-restrictive smoothness 
assumptions for the case of i.i.d. random variables. In the present paper we extend his main result to 
the case of serial dependent random variables by means of the powerful and elegant functional delta 
method. Moreover, we utilize the functional delta method in order to obtain conditional consistency of 
certain bootstrap procedures. Finally, we extend the results to the more general sequential empirical 
copula process under serial dependence. 
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1 Introduction 



x 



Let F be a d-variate continuous distribution functio n, d > 2, wit h marginal distribution functions F p for 
p E {1, . . . , d}. By the famous Theorem of Sklar [see ISklarl ([1959T )] we can decompose F as follows 



F(x) = C(F 1 (x 1 ),...,F d (x d )), x= ( Xl ,...,x d ) e 



1.1 



where C is the unique copula associated to F. By definition, C is a <i-variate distribution function on 
the unit cube [0, l] d whose univariate marginals are standard uniform distributions on the interval [0, 1]. 
Equation (|1.1|) is usually interpreted in the way that the copula C completely characterizes the information 
about the stochastic dependence c ontained in F . For an extensive exposition on the theory of copulas we 
refer the reader to the monograph iNelsenl (|2006l ). 
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A fundamental tool for statistical inference on copulas is the so-called empirical copula. If (Xj)j g z forms 
a strictly stationary sequence of d-variate random vectors with copula C and if F n denotes the empirical 
distribution function of Xi, . . . , X n , then the empirical copula is defined as 



C n (u) = F n (F~ 1 (u 1 ),...,F~ d (u d )), u€ [0,1]° 



where F np and F np denote the p-th. marginal empirical cdf and its generalized inverse, respectively, for 
p = 1, . . . , d. The asymptotic behavior of the corresponding empirical copula process 



C„(u) = V^(C n (u) - C(u)) 



was studied for i.i.d. samp l es in 



Gaenssler 



and 



Stutel (119871 ): 



Fermanian et al 



dgOOJ ); iTsukaharal ((2003 ): 



van der Vaart and Wellnerl (|2007l ): ISegera (|201ll ). among others. Except for the results by the last named 
author, all previous deri vations are b ased on the assumption that C has continuous partial derivatives 
on the closed unit cube. ISegersI (|201ll ) pointed out that there are hardly any popular copulas satisfying 



this condition and showed that the results actually do hold under much less restrictive assumptions which 
are fulfilled for many co pula models t hat are applied in practice. As we will demonstrate in Section [21 
the assumption made by Segers (|2011 ) is in fact sufficient for weak convergence of the empirical copula 
process under various kinds of weak dependence. 

In the case of weakly dependent str ictly stationary s ampl es Xi , . . . , X n only few results are available to 
date. To the best of our knowledge, iDoukhan et al.l (|2009i ) is the only reference where weak convergence 



of the empirical copula process is established for a sample of //-dependent random variables; again under 
the restrictive smoothness assumption of continuous partial derivatives on the whole unit cube disc ussed 
above. It is on e of th e main purposes of the present paper to extend the results of ISegersI (|201ll ) and 



Doukhan et al 



(2009) to a weak convergence result for C n under both a non-restrictive smoothness 
assumption and under various weak dependence concepts. At the same time, we provide simple conditions 
that allow for future generalizations of results on the empirical copula pr o cess. Our method of proof 
is based on the functional delta method [see also 



Fermanian et al 



(2004); 



Van der Vaart and Wellner 



(1996)], which implies that the weak convergence of the classical empirical process combined with a 
natural property of the limiting distribution and the mild smoothness assumption of ISegera (120111) is 
alread y sufficient for the weak convergence of C n . As a by-product we rediscover the findings from lSegers 
(2011) in a way that circumvents sophisticated probabilistic arguments. Moreover, the functional delta 



method allows to prove consistency of bootstrap procedures for the empirical copula process that are 
constructed from consistent bootstrap approximations of the classical empiri cal process. 



The second purpose of the present paper is an extension of the results of iRuschendorfl (|1976l ) for the 
sequential empirical process. More precisely, we consider weak convergence of the following two processes 



[sn] 

-=£(l{X 4 <F-(u)} 

^ n i=i 



[sn] 



C(u)) 
<F-(u)}-C n (u)) 
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where s E [0, 1] and F n (u) = (F n x(ux), . . . , F nd {ud)). The first process is the one considered in lRiischendorf 
( 19761 ) [except for a small error term of order n _1 arising from the fact that F~ (F np (u p )) is not exactly 



equal to u p ] and it contains the usual empirical copula process as a special case, since or(i , u) = C n (u). 
As in the case of the usual empirical cop ula process C n described above, we both weaken the smoothness 



assumptions made in 



Riischendorf 



(|1976l ) and we allow for various weak dependence concepts in order to 



i —n 

Wied et al. 


(2011 


); 


Remillardl ( 


2010) or 


Ruppert 


(2011) 



obtain the weak convergence of or. 



constant dependence over time, see e.g. 
for the last named reference, the asymptotics of are derived from the ones of Cjf • However, the 
limiting process of C£ indicates that those assumptions should actually not be necessary. The latter 
author shows for the case of strictly stationary samples that C+ indeed converges weakly without any 
smoothness assumption on C at all. In Section [3] we treat this problem by the functional de lta method, 
which easil y allows fo r more general weak dependence concepts. Note, that the methods in 



Wied ct al. 



(2011 



and 



Remillardl (|2010l ) are greatly improved in their applicability by our results. 



The paper is organized as follows. In Section [2] we introduce the basic notation for the remaining part 
of the paper. Moreover, the usual empirical copula process, its weak convergence and a (blockwise) 
bootstrap approximation under serial dependence is investigated. Some simple examples of applications 
of our results are given for illustration. In Section Owe derive similar results for the sequential empirical 
copula processes. 



2 The empirical copula process 

2.1 Weak convergence of the empirical copula process 

Let (Xj)j g z with Xi = (Xn, . . . , Xid) be a strictly stationary sequence of <i-variate random vectors. 
Throughout the paper, the joint distribution function F of Xj is assumed to have continuous marginal 
distributions Fx, ■ ■ ■ and its copula is denoted by C. Further, let Ui P = F p (Xi p ) be the probability 
integral transform of the components of X, and set U} = (Un, . . . , Uid). Note that Uj is distributed 
according to C . Consider the empirical distribution functions 

n n 

F n (x) = n" 1 ^I{X, < x}, G n (u) = n- 1 ^I{Ui < u} 

i=l i=l 

for x G M. d , u G [0, l] d . For some distribution function H on the real line let 

H-(p) :- 



mf{x E M | H{x) > p}, 0<p<l 
sup{x E R | H(x) = 0}, p = 



denote the (left-continuous) generalized inverse function of H. The object of interest is the empirical 
copula defined as C n (u) = F n (F~ 1 (ui) , . . . , F~ d (ud)) , where F np denotes the p-th marginal empirical 
cdf, p = 1, . . . , d. For convenience, we will abbreviate the notation for the empirical copula by C n (u) = 
F n (F~(u)), with F~(u) = . . . , F~ d (ud)). Note that G n can not be computed from the data unless 
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the margina l distr ibutions Fi, . . . , Fd are known. However, it can be shown that C n = G n (G~ 1 , . . . , G~ d ) 
[cf. ISegera (|201ll )]. which is the reason why we can prove the subsequent theoretical results by just 
considering the sequence (Uj)j e z and its empirical distribution function G n . 

In order to derive the asymptotics of the corresponding empirical copula process C„ = y/n(C n — C) we 
make the following assumption on the usual (cf-variate) empirical process a n = y/n{G n — C) based on 
Ui,...,U n . 

Condition 2.1. For the strictly stationary sequence (Xj)j g z it holds that a n = y/n(G n — C) weakly 
converges in ^°°([0, l] d ) to a tight, centered Gaussian field Be concentrated on Do> where 

Do = | a £ C[0, l] d | a(l, . . . , 1) =0 and a(x) = if some of the components of x are equal to o| . 

Here, and throughout the rest of this paper £°°{[0,l] d ) denotes the space of all uniformly bounded functions 
on the unit cube [0, l] d , equipped with supremum norm || • ||oo* Note that Condition [2.11 is a standard 
result of empirical process theory for weakly dependent data, where the limiting process Be is a tight, 
centered Gaussian process on [0, l] d with covariance 

r(u,v) = Cov(B c (u),B c (v)) = ^Cov(I{U < u},I{U, < v}). 

jez 

Sufficient conditions that entail Condition 12.11 have been proposed in terms of various weak dependence 
concepts: 



for q- mixing (str ongly mixing) conditions, see Theorem 2 in lPhilipp and Pinzurl (|1980l ) or Theorem 



7.3 in 



Rio (2000) for a more refined version. 



for /?-m ixing (absolutely regular) conditions, see e.g. Chapter 11.6 in lKosorokl (|2008l ) or 
([l995h . 



Doukhan et al. 



for ^-dependence, see Theorem 1 in 



Doukhan et al 



(2009). 



for multiple mixing properties see lDehling and Durieul (|201ll ) or lDurieu and Tuschd (|201ll ). 



for a similar result under long range dependence [at a different rate, though] see iMarinuccil (|2005l ) 
and Remark 12.61 a) below. 



It is the main purpose of this paper to prove weak convergence of C n under the general Condition 12.11 
incorporating all the aforementioned citations. 

Example 2.2. (ARMA(p, g)-processes) The Revalued process (X t ) tg ^ is said to be an ARMA(p, q)- 
process if it is stationary and if it satisfies the equation 

v Q 



i=0 



j=0 



for every t G Z, where B(i) G R dxd with B(0) = Id, A(k) € R dxr (r > 1) and where Z t G W are i.i.d. 
centered random vectors with covariance matrix E. For z G C let P(z) = Yli=o B{i)z % . It was shown in 
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Mokkademl (|l988l ) that if the absolute values of the zeros of the polyomial det P{z) are strictly greater 
than 1 and if the probability law of Z t is absolutely continuous with respect to the Lebesgue m easure on 
, then Xt is /3-mixing with f3(k) = 0(e~ ak ) for some a > 0. For instance by the results in 



Kosorok 



(|2008l ). this is sufficient for Condition 12,11 to hold, with much to spare. 



For various further examples of stochastic processes (e.g. GARCH processes, and various extensions 
thereof) satisfying any of the weak dependence concepts noted above we refer the reader to the corre- 
sponding citations. 

In order to accomplish the step from the known weak convergence of a n to the weak convergence of C n 
some smooth ness assumptions on C are needed [even in the simple point- wise sense]. As pointed out by 
Segersl (|201ll ) the following smoothness condition on C is nonrestrictive in the sense that it assures that 
the candidate limiting process of C n has continuous trajectories. 

Condition 2.3. For p = 1 . . . , d the partial derivatives d p C(u) exist and are continuous on U p = {u € 
[0,l] d \u p G (0,1)}. 



Now consider the mapping 



£°°[0,l] c 



H i-> H(H 



l ' ■ 



H(H- 



where B<j> denotes the set of all distribution functions H on [0, l] d whose marginal cdfs H p satisfy H p (0) = 
0. Observe that 

C n = V^(C n -C) = v^(G n (G-, . . . , G7 nd ) -C) = V^mGn) - $(C)). 

Therefore, in order to derive the weak convergence of C n it suffices to establish Hadamard-differentiability 
of $ tangentially to suitable subspaces and to invoke the fun ctional delta method. This Hadamard- 

()2004l ) und er the restrictive assumption that C 



Fermanian et al. 



differentiability has been investigated in 

has continuous partial derivatives on the whole unit cube [0, l] d [cf. IVan der Vaart and Wellnerl (| 19961 ). 
page 389] . The following Theorem relaxes this assumption to Condition 12.31 its proof is given at the end 
of this section. 

Theorem 2.4. Suppose Condition \2.3\ holds. Then <I> is Hadamard-differentiable at C tangentially to 
Bo- Its derivative at C in a £ Bo * s given by 

d 

P =l 

where d p C is defined as on the set {u € [0, l] d \ u p G {0, 1}} and where vS-^ = (1, . . . , 1, u p , 1, . . . , 1) 
denotes the vector, where all entries of u except the p-th are replaced by 1. 



The functional delta method, see e.g. Theorem 3 .9.4 in Van der Vaart and Wellnerl (|1996 ), immediately 



yields the following result. It has been stated in ISegerd (|201ll ) for the i.i.d. case and in 



Doukhan et al. 
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(2009) for the serial dependent case, where the latter authors impose restrictive smoothness conditions 
on C [and investigate the Skorohod-space D([0, l] d ), which implies weak convergence with respect to the 
sup-norm if the limiting process has a version that has continuous trajectories almost surely]. 

Corollary 2.5. Suppose Conditions \2.1\ and \2.3\ hold. Then the empirical copula process C n = y/n(C n — 
C) weakly converges towards a Gaussian field Gc, 



which may be expressed as 



C n = Vn~(C n -C) ^ G c in P°([0,l] d ), 



G c (u) = B c (u) - ]T<9 p C(u) Bc(u (p) ; 
P =i 



(2.1) 



Remark 2.6. 

a) More generally, the functional delta method allows the derivation of weak convergence results under 
different rates of convergence of the em pirical process. F or instance, in the case of long range dependent 



stationary sequences, it was shown by 
converges in -D([0, l] 2 ) to some limit process 



Marinuccil (J2005) in the bivariate case that r n (G n — C) weakly 
£ Bo, where the rate of convergence l/r n — > satisfies 



o jn 1 / 2 ). 



f fut ure research reveals more general results on the usual empirical process than the one 



by iMarinuccil (|2005l ) our method of proof easily extends these results to the empirical copula process. 



b) Note that a direct extension of the method of proof employed by ISeeerd (|201ll ) to the case of weakly 
dependent stationary sequences would need additional probabilistic arguments that differ depending on 
the underlying data structure. On the other hand, the functional delta method allows to separate the 
probabilistic part (weak convergence of the empirical distribution function and properties of the limiting 
process) and the analytical part (compact differentiability) of the problem hereby considerably simplifying 
the proof of the weak convergence of C n . 



The following example is a simple, but typical application of Corollary 

Example 2.7. (Asymptotics of Spearman's rho) 

The multivariate population version of Spearman's rho is defined as 



P 



d+1 



2 d -(d + 



i)( 2i L c<u)< H' 



sec 



Schmid and Schmidtl (|2007l ). and may be estimated by the simple plug- in nonparametric estimator 

d + 1 /„ 



Pn 



C n (u) du — 1 



2^- (d+1) \ j md 
Observing continuity of the mapping C \-+ p(C) we immediately obtain the result that 

d+1 



2 d -(d + 1) 



[0,l] d J [0,1] 



E [G c (u)Gc(v)] dudv. 
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2.2 Bootstrapping the empirical copula process 

In the present section we apply the functional delta method in order to obtain consistent bootstrap 
procedures for the empirical copula process under a general assumption on a given bootstrap procedure 
for the usual d-variate empirical process. As an illustration of our results, we provide consistency of a 
simple (blockwise) bootstrap approximation of the limiting field Gc in the setting of a-mixing data. More 
precisely, let Mi, . . . , M n be non-negative, real- valued random variables with Ya=i M% = n and define 

n n 

F n , 6 (x) = n- 1 Mi I{X, < x}, G n>6 (u) = n' 1 M t I{Uj < u}. 

i=l i=l 

For instance, if (Mi, . . . , M n ) follows a multinomial distribution with parameters n,n , . . . ,n , we 
obtain the usual resampling bootstrap of an i.i.d. sample, i.e. the bootstrap based on sampling from 
the data with replacement. On the other hand, the construction above also allows to handle bootstrap 
procedures such as the blockwise bootstrap considered in Example 12.101 below. 

In order to construct a consistent bootstrap procedure for the empirical copula process, we need the 
following condition. 



Condition 2.8. Additionally to the assumptions made in Condition 12. 1\ suppose that a„& = y/n(G n b ■ 
G n ) weakly converges in £°°([0, 1] ) to Be conditional on the data in probability, notationally a n ^ ~» B( 



Weak convergence conditional on the data in probability ( -convergence) is understood in the Hoffmann- 



J0rgensen sense as defined in Kosorok ( 2008 ^. that is a n ^ ~~* a if and only if 



M 

(i) sup /gBLl |Ejv//(a n)b ) - E/(a)| ->• in outer probability, 

(ii) E M f(a n>b )* - M M f(a n)b )* 4 for all / € BL l5 

where BLi denotes the set of all Lipschitz-continuous functions / : £°°([0, l] d ) — > R that are uniformly 
bounded by 1 and have Lipschitz constants bounded by 1, and where the asterisks in (ii) denote measurable 
majorants (and minorants, respectively) with respect to the joint data (Xi, . . . , X n , Mi, . . . , M n ). We give 
three examples under which Condition 12.81 is met. 

Example 2.9. (Bootstrap in the i.i.d. case) 

a) Suppose Xi, . . . , X n are i.i.d. Moreover, assume that £i, . . . , are i.i.d., non-negative, real- valued ran- 
dom variables that are independent of the Xj and satisfy the conditions E£i = l,E£f = 1 and ||£i||2,i = 
J °° \/P(|^i| > x) dx < oo. D efine the rando m variables Mi,...,M n as Mi := ^/(^" =1 ^),z = 



1, . . . , n. Then Theorem 2.6 in 



Kosorokl (j2008l ) implies that Condition 12.81 holds. 



b) Assume that Xj, . . . , X n are i.i.d. and that the vector (Mi, . . . , M n ) f o llows a multinomial distribution 
with parameters n, n _1 , . . . ,n _1 . Again, by Theorem 2.6 in iKosorok (|200al ) Condition!^ holds. 
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Example 2.10. (Bootstrap in the strongly mixing case) 

Suppose (Xj)igz is a strictly stationary, strongly mixing sequence whose a-mixing co e fficien ts satisfy 
the summability condition £°^ (i + l) 16(d+1) a 1/2 (i) < oo as considered in iBiihhnannl (|1993l ). In this 
setting, the usual bootstrap based on sampling with replacement from the data (see Example 12.91 b)) 
is not able to ca pture the serial dependence of the sequence. The blockwise bootstrap considered in 
Buhlmannl (119931 ) provides one possible solution. It is based on the following idea. Instead of sampling 
(with replacement) single data points, we consider n — I + 1 blocks of length I which have the form 
(Xi, . . . , X;), (X2, • • • , X^ +1 ), . . . (X n _; + i, . . . , X n ). The bootstrap data are formed by sampling with re- 
placement k = n/l of those blocks. This works for k being an integer, if that is not the case we simply 
make the last selected block smaller. If I = l n is chosen appropriately, this procedure is able to capture the 
serial dependence structure in the data. More formally, we define the following bootstrapped empirical 
distribution functions [for the sake of a clear exposition we only consider the case k = n/l € N] 

k 



k Si+l 

F nib (*) = k- 1 Y j r 1 j2 i{Xj <x}, 

i=l j=Si+l 



G n , b (u) = k 



i=i 



Si+l 

E HU^u}, 

j=Si+l 



where n = kl,l = l n = o(n), l n — >• 00 and (Si)i^ are i.i.d. uniformly distributed on the set {0, . 



,n 



Note that in the case n = kl the quantity F n &(x) can also be represented as F n b (x) 



n 



x} where Mj := #{i = 1, . . . , k \ j — I < S{ < j — 1}. Thus we are in the setting describe d at the be 



gmnm 



Biihlmann 



(|1993l l 



of this section. Given that l n — > 00, l n = 0(n l / 2 ~ £ ) for some e > 0, Corollary 3.1 in 
states that a n ^ ~^ Be almost surely in (D([0, l] d ), \\ ■ \\oo)j where the latter space is equipped with 
the cr-field induced by the open balls. Note that a n ^ is measurable with respect to the ball-u-field and 
the theory for weak convergence on metric spaces spaces equipped with non-Borel sigma fields [see e.g. 
Pollard! ()1984l ). Chapter IV] applies. Thus the delta method for weak convergence in the Hoffmann- 
J0rgensen sense is not directly applicable. Howev er, a closer look at Theorem 1.7.2 and the proof of 



Van der Vaart and Wellnerl (jl996l ) reveals that we can enlarge both the space -D([0, l] rf ) 

V^(Gn,b-G n ) -w M c inf°([0,l] d ) 



Theorem 2.9.6 in 

to £°°([0, l] d ) and the cr-field to the Borel cr-field and obtain a n ^ 
in the Hoffmann-j0rgensen sense. 



Now, defining 

C n , b (u) = F n , b (F- b (u)) = G njb (G- b (u)) 

as the b o otstra p empirical copula, the functional delta method for the bootstrap, see Theorem 12.1 in 
Kosorokl (|2008l ). yields the following result. 



Corollary 2.11. 



Under condition 1 2. 3\ and 

\At(c„ jb 



it holds 



C n )^G c in r°([0,l] d ). 

M 



This result has also been stated in 
a) and in 



Fermanian et al. 



Biicher and Dettd (|2010l ) under the special setting of Example 12.9 



(2004) under the setting of Example 12.91 b), where both citations use the 
restrictive smoothness assumption that the copula has continuous partial derivatives on the whole unit 
cube. 
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Remark 2.12. In the i.i.d. case the literature provides an alternative multiplier bootstrap pr o cedur e 
based o n the estimation of the partial derivatives of the copula, see e.g. iRemillard an d Scaillet (|2009l ): 
Segers ( 2011 ). Consistency is usually stated in a weaker unconditional sense [see Bucherl ( 2011 ) for an 
extension to the conditional Hoffmann- J0rgensen sense]. The investigation of this bootstrap procedure 
in the serial dependent case is beyond the scope of the pre s ent pa per since it cannot be handled by the 
functional delta method. We refer the reader to iRuppertl (|201ll ) for the consideration of the strongly 



mixing case. 



Finally, we consider one of many possible applications of Corollary 12.111 

Example 2.13. (Testing for symmetry of a copula) 

Suppose one is interested in a test for the (bivariate) hypothesis Hq : C(u, v) = C(v, u) for all u,v G [0, 1]. 
Recently, 
process 



Genest ct al 



(|201ll ) proposed to base a test for T-Lq by considering functionals of the following 



B n (u,v) = \/n(C n (u,v) - C n (v,u)) = J C n (u,v) - C n (v,u) 

which weakly converges to O c (u,v) = G c (u,v) - G c (v,u) in £°°(\0 l} 2 ) by the continuous mapping 
theorem. By the bootstrap continuous mapping theorem, see 



Kosorok 



(|2008h . we obtain 



lt b(u,v) = Vn(C njb (ii, v) - C n! b(v,u)) ~> B c (u,v) 



independently of the hypothesis. One test statistic investigated by 



Genest et al 



(|201ll ) is defined as 

T n = 1 1 D n | |oo, and we reject Ho if T n is larger than the 1 — a quantile of the law of Tq = || Del loo; which 
depends on the unknown copula in a complicated manner. Approximate quantiles can be obtained by a 
bootstrap sample of T n ^ = ||B n b|| 00 . 



2.3 Proof of Theorem EH 

Let E denote the set of distribution functions F on [0, 1] with .F(O) = and define E~ as the set of all 
generalized inverse functions F~ with FgE. Now decompose <3? = $3 o $ 2 where 



$9 



H ^ {H,H l ,...,H d ) 

D)$ 2 B$ 3 = B$ x (E-) c 

(H,F 1 ,...,F d ) ' y (H,Fi,...F-) 

r°([o,i] d ) 



"*3 



(H,G 1 ,...,G d ) ^ Ho(G 1 



,G d ). 



(2.2) 



The first mapping <I>i is Hadamard-differentiable at C with derivative = ^1 since it is linear and 
continuous. 

Considering the second mapping let U G E be the identity, then A : F 1— > F~ with F £ E is Hadamard- 
differentiable at U tangentially to the set E = {F £ C[0, 1] : F(0) = F(l) = 0} with derivative 
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M v (h) = -h. To see this let t n -> 0, h n G ^°°([0, 1]) with li n -> /i e E and *7 + t„/i„ G E. We have to 
show that 



sup 

xG[0,l] 



A(U + t n h n ) - A(U) 



A' v (h) 



(x) = sup 

££[0,1] 



(U + t n h n ) (x)-x 



+ h(x) 



0. 



(2.3) 



We deal with the case x = first and show that 

A, = (U + i,A»)~(0) = sup{x G [0, 1] : x + t„/i„(x) = 0} = o(t n ). 

Suppose ^4 n is non-zero, eventually. Choose a sequence x n G [A n /2,A n ] satisfying x n + t n h n (x n ) = 0, i.e. 
h n (x n ) = —Xn/tn- Since h n is bounded we obtain x n = 0(t n ) = o(l). By uniform convergence of h n and 
continuity of h this implies —x n /t n = h n (x n ) — > h(0) = 0, that is x n = o(t n ) and therefore A n = o(t n ). 
Now consider the case x > 0. Since U + i n ^n is a distribution function with ([/ + t n h n )(0) = it follows 
that £ n (x) = (U + t n h n )-(x) G (0, 1] for all x G (0, 1]. Set e(x) = A > 0, then 

(U + t n h n )(£ n (x) - e{x)) <x<(U + t n h n )(£ n (x)) 

for all x G (0, 1] by the definition of the generalized inverse function. This implies 

-t n h n {tn{x)) < 5,n{x) ~X< -t n h n (^ n {x) - e{x)) + t\. 

Since h n is uniformly bounded we obtain £ n (x) — > x uniformly in x G (0,1]. Divide the last equation 
by t n and use uniform convergence of h n and continuity of h to conclude that (Cn(x) — x)/t n — > —h(x) 
uniformly in x G (0, 1]. Together with the = this yields (j2.3|) . 

It now easily follows that $2 is Hadamard-differentiable at (C, U,...,U) tangentially to Bo x Eq with 
derivative 



2,(C,U,...,U) 



(a, fx, ...,fd) = («, -fx, ■■■ , ~fd)- 



Last but not least consider $3. We assert that $3 is Hadamard-differentiable at (C, U, . . . ,U) tangentially 
to Do x Eq with derivative 

d 



ft 



3,(C,u,...,U)( a i9i, ■ ■ ■ ,9d)(*) = "( x ) +J2d p C(-x.)g p (x p ) 



p=X 

To see this let t n and {a n , g nl , . . . , g nd ) G r°([0,l] d ) x (f>°([0, l])) d with ( 

<^n > 9nX 5 • • • 5 9nd ) ' 

(a,gx, • • • , gd) G Bo x Eq such that (C + t n a n , U + i„5„i, . . . , U + t n g nd ) G P$ x (E~) d . Now decom- 
pose 



where 



t n 1 { $ 3(C + t n a n , U + i n # nl , ...,U + t n g nd ) - $ 3 (C, 17, . . . , U)} = L nl + L n2 



L m = t^iCo (U + t n g n x, ...,U + t n g nd ) - C} 
L n 2 = a n o (U + t n g n i, ...,U + t n g nd ) 
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and consider both terms separately. Exploiting the facts that a n and (g n i, . . . ,g n d) converge uniformly 
and that a is uniformly continuous one can conclude that \\L n 2 — a||oo = o(l). Concerning L n \ we have 
to deal with different cases regarding the number of coordinates on the boundary of [0, l] d . If x € (0, l) d 
a Taylor expansion of C at x yields 

d 

L„i(x) = ^2 dpC(x.)g np (x p ) + r n (x), 
p=i 

where the error term r n can be written as 

d 

r «( x ) = ^2 {9 P C(y) - d p C(x))g np (xp) 
P =i 

with some intermediate point y = y(n) between x and (U + t n g n i, . . . , U + t n g nd ){x). The main term 
uniformly converges to J2p=i 9pC{x.)g p (x p ) as asserted, hence it remains to show that the error term 
converges to uniformly in x. To see this, let e > 0. Using uniform convergence of g np , uniform 
continuity of g p and the fact that g p (0) = g p (l) = one can conclude that there exists 5 > such that 
\gnp(x p )\ < e/2 for all x p < 8 and x p > 1 — 5 and all p = 1, . . . , d. Since partial derivatives of copulas are 
bounded by 1 the p-th. summand of the error term is uniformly bounded by s for x p < 5 and x p > 1 — 5. On 
the set [0, l] d \ {x p < S or x p > 1 — 6} the partial derivative d p C is uniformly continuous which yields the 
desired convergence under consideration of y(n) — > x and boundedness of g p . Since e > was arbitrary 
the case x G (0, l) d is finished. 

We introduce the notation X( p ) = (x±, . . . , x p -±, l,x p+ \, . . . , x^) and now consider the case x = X( p ) with 
x q G (0, 1) for q p. Decompose L n i(x) = L^(x) + L^(x) where 

L nl ( X ) = *n 1 { C '(( a; l + *nffnl(a;i), ■ • • ,X d + t n g nd (x d ))( p) ) - C(x (p) )} 
(2) —1 

L ni ( x ) = l n {C{x\ + t n g nl (xi), ...,x d + t n g nd (x d )) - C((xi + t n g n i(xi), ...,X d + inffnd^d))^))}- 

By the same arguments as in the previous case a Taylor expansion in X( p ) and a discussion of the remainder 
yields 

d 

L£(x) = J> g C(x) 5g (z g ) + = J2d q C^)g q (x q ) + o(l) 

(2) 

uniformly in x (note that g p (l) = 0). Lipschitz-continuity of C implies \L nl (x)| < |^ n p(l)| — > |^ p (l)| = 
and the second case is finished. If more than one coordinate is 1 while the others stay in (0, 1) the 
argumentation is analogous. 

If one (or more) of the coordinates of x equals 0, say x p = use C(x) = and Lipschitz-continuity of C 
to estimate 

|im(x)| < \g np (,0)\ ->0. 

If Xp = we have d q C(x) = for q / p and g p (x p ) = g p (0) = 0, which entails = Ylp=i d p C(x)g p (x p ) 
and thus the assertion. 
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The remaining case x = (1, . . . , 1) again follows by Lipschitz-continuity, the details are omitted. To 
conclude, $3 is Hadamard-differentiable as asserted. 

Now apply the chain rule, Lemma 3.9.3 in IVan der Vaart and Wellnerl (| 19961 ). to 3> = $3 o $ 2 o $1 to 



conclude the assertion of the Lemma. 



□ 



3 The sequential empirical copula process 

Under the setting of Section [2] consider the following partial sum 

[sn\ 

G#( S ,u) = n- 1 ^I{U i <u} 

i=l 

for s G [0, 1] and u E [0, l] d . The process 

[sn\ 

a#(s, u) = n~V2 £ ( I{Ui < u} - C(u)) = ^ (g#(s, u) - C#(s, u)) + Ofr- 1 / 2 ), 
i=i 



Van der Vaart and Wellner 



where C*(s, u) = sC(u), is called sequential empirical process, see e.g. Chapter 2.12 in[ 
(1996). The remainder term 0(n -1 / 2 ), holds uniformly in s,u and over all copulas C. Estimating the 
marginal distribution functions in a„ by their empirical counterparts we arrive at the so-called sequential 
empirical copula processes defined as 

#( S ,U) = -=Y. ^ F n(")} ~ C ( U )) = ^ ~ ^C(U)J 

i=l 

= v^(c#(s,u) -C # (s,u)) +0(n" 1 / 2 ) 

1 ^ / II \ 

C+( S ,u) = -= (l{Xi < F~(u)} - C n (n)) = [C*(s,u) - i^C#(l,u) j 

i=i ^ ' 

= (C*{s, u) - a(7 n (u)) + ¥ {n-^ 2 ) 

where F~(u) = (i^fai), • • • , F~ d (u d )) and C#(s, u) = n" 1 I{X< < F~(u)}. For s = 1 we retrieve 

the usual empirical copula proce s s from Section [2j i.e. Cn(l, •) = C n . The process Cn was introduced 
and investigated in iRiischendora ()1976l ) . In order to state his general result we will make the following 
assumption on the weak convergence of the sequential empirical process at, which is similar to Condition 
2Tlin Section El 



Condition 3.1. For the strictly stationary sequence (Xj)j £ 2 it holds that at = y/n(Gt — C*) weakly 
converges in £°°([0, l] d+1 ) to a tight, centered Gaussian field concentrated on Dq , where 



4 



I a G C([0, l] d+1 ) I a(0, •) = and a(s, •) G B for all s G (0, 1]} . 
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As for the usual empirical process, the literature provides several sufficient conditions for the weak con- 
vergence postulated in Condition 13.11 For instance, in the case of i.i.d. random vectors we have weak 
convergence to the Kiefer-Miiller process, which has mean zero and covariance function 

Cov(B#(s,u),B#(t,v)) = (sAt) (C(uAv) - C(u)C(v)) , 



see e.g. Chapter 2.12 in | Van der Vaart and Wellnerl ()1996l ). For a similar result under strongly mixing 
conditions see e.g. Iphilipp and Pinzurl Jl98fT i" [which again incorporates the ARMA(p g)-process discussed 



in Example 12.21 ] 

For technical reasons we consider the asymptotic behavior of C+ first. For G £°°([0, l] d+1 ) let -ff(x) = 
i?*(l,x) and consider the mapping 

-> r°([0,l] d+1 ) 
H* i — y ^(F # )(s,x) = H#(s,H-(x)) - sH(H~(x)) 

where B# consists of all functions H* £ £°°([0, l] d+1 ) such that H(-) = H*(l,-) £ B$. Recall from 
Section [2] that C n (u) = G n {G~(u)). Similarly, it holds that Cjf (s,u) = Gn(s,G~(u)) and thus 




s,u) = V^[G*(s,G-(u)) - sG n (G~(u))) +O p ( 



ii 



(*(G#) - *(C # )) (s, u) + P (n" 



where the last equality follows from the fact that ^(C*)(s,u) = C*(s,u) — sC*(l,u) = 0. Similarly, 
define 



B* -»■ r°([0,l] d+1 ) 

F# i — ^ r(tf#)(s,x) = H*(s,H~(x)) 



r : 



and note that C« = \/n ^T(G>f ) — r(C*)j + Op(n 1 / 2 ). Mimicking the argumentation in Section [2j in 

order to derive the asymptotics of and C|f it remains to show Hadamard-differentiability of \& and 
r at [tangentially to suitable subspaces, if necessary]. This is done in the following Theorem. Note, 
that for the first part we do not need any assumptions on the smoothness of C, the result holds for every 
copula C. 

Theorem 3.2. a) Let C be an arbitrary copula. Then the mapping ^ is Hadamard-differentiable at 
6 B\j,, C*(s,x) = sC(s,x) with derivative 

^ # (a # )(s,x) = a # (s,x) - sa # (l,x). 

b) If moreover the copula satisfies condition \2.3l then the mapping F is Hadamard-differentiable at C# 
tangentially to B* . Its derivative at in oft 6 B* is given by 

d 

T' c# (a*)(s, x) = a#(s, x) - s £ d p C(x) a*(l, x«). 

p=i 
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Proof. We begin with the proof of a). Let t n -> and a$ G £°°([0, l] d+1 ) with af -> a # G r°([0, 
such that C # + f n af G D*. Then 

= ^{(C* + t n a*)(s, (C + t„a n )-(u)) - s{C + t„a n ))((C + t n a n )-(u))} 
= a*(s, (C + t n a n )~(u)) - sa n {(C + t n a n )~(u)). 

By the proof of Theorem 12.41 we can conclude that (C + t n a n )~ uniformly converges to the identity on 
[0, l] d , see equation (|2.3p . whose proof did not need Condition 12.11 Together with uniform convergence of 
0$ and uniform continuity of aft we obtain assertion a). 

For the proof of b) note that r(a*)(s,x) = ^(oft)(s,x) + s$(a)(x) with $ and $ from Theorem 13.21 
and Theorem 12.41 respectively. Some thoughts reveal that the Hadamard-differentiability of ^ and $ 
transfers to T, such that r^ # (a*)(s,x) = $^, # (a*)(s,x) + s < I )/ C (a)(x). The details are omitted for the 
sake of brevity. □ 

The following main result of this section is a consequence of Theorem 13.21 and the functional delta method. 
Again note, that we do not need any smoothness assumptions on C for the first part of the result. 

Corollary 3.3. a) Suppose Condition \3. 1\ holds. Then we have, for any copula C, that 

C+( S ,u) = v^(c#(s,u)- S C n (u)) +0 P (n- 1 / 2 ) ~» G+ in l°°([Q, l] d+1 ), 

where the limiting Gaussian field may be expressed as Gq(s,u) = B*(s,u) — sB^(l,u). 
b) If moreover the copula satisfies Condition \2.3\ then we have 



C#( S ,u) = v ^(c#(s,u)- S C(u))+0 P (n~ 1 / 2 ) ^ G* in Z°°([0, l] d+1 ), 
where G*{s,u) = M*(s,u) - d p C(u) B#(1, U W). 

Corollary 13.31 has several applications for the investigation of multivariate tim e series, wh e n one is inter- 

) test for 



Wied et al 



(2011 



ested in tests for structural breaks in the dependence structure. For instance, 
a constant Spearman's rho over time of some strictly stationary time series, while iRemillardl (|2010l ) tests 
for a constant copula. The results of both papers are improved by Corollary 13.31 a) , since we do not make 
any smoothness assumption on C at all. 
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